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j> to investigate supersymmetry on the Lattice. 

• ^ 

X 



* antipin@cp3-origms.net 



moiaza@cp3-origins.net 


pica@cp3-origins.net 




sannino@cp3-origins 


net 



1 



I. INTRODUCTION 



Gauge theories constitute the building blocks of our present understanding of natural 
phenomena. The Standard Model (SM) of high energy particle interactions is, in fact, 
entirely based on a semi-simple gauge group SU{3) x SU{2) x U{1). The SM accounts 
for roughly four percent of the known universe. However, several puzzles remain still 
unexplained. For example, why do we observe, at least, three generations of elementary 
particles. Besides, the remaining 94% of the universe, of which 22% dark matter, and 
72% dark energy, remains largely unknown. It is therefore natural to expect that new 
gauge theories, or extensions of the present ones play a fundamental role in explaining 
the unknown side of the universe. The space of four dimensional gauge theories at our 
disposal, without theoretical prejudice, is very large and moreover a large fraction of it is 
still terra incognita due to our limited methods to tackle nonperturbative dynamics. 

A fascinating possibility is that different gauge theories display the same physical 
dynamics. This is possible since the physical quantities are determined by the gauge 
singlet operators of a generic gauge theory. Therefore it is perfectly legal to imagine to 
have two different gauge theories leading to the same physics, at least, in certain regimes. 
Quantum chromodynamics (QCD) itself is a famous example. In fact, at low energy 
the spectrum can be represented by hadrons interacting via Yukawa interactions and 
its electric dual description is in terms of the SU{3) gauge theory featuring quarks and 
gluons. 

It is therefore natural to imagine that generic asymptotically free gauge theories have 
magnetic duals. In fact, in the late nineties, in a series of ground breaking papers Seiberg 
im 121 provided strong support for the existence of a consistent picture of such a dual- 
ity within a supersymmetric framework. Supersymmetry is, however, quite special and 
the existence of such a duality does not automatically imply the existence of nonsuper- 
symmetric duals. One of the most relevant results put forward by Seiberg has been the 
identification of the boundary of the conformal window for supersymmetric QCD as a 
function of the number of flavors and colors. Recently several analytic predictions have 
been provided for the conformal window of nonsupersymmetric gauge theories using 
different approaches Il3]-|3T1. We initiated in ||32| the exploration of the possible existence 
of a QCD nonsupersymmetric gauge dual providing a consistent picture of the phase di- 
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agram as a function of number of colors and flavors. Arguably the existence of a possible 
dual of a generic nonsupersymmetric asymptotically free gauge theory able to repro- 
duce its infrared dynamics must match the 't Hooft anomaly conditions ||33|. We have 
exhibited several solutions of these conditions for QCD in ||32l. An earlier exploration 
already appeared in ||34|. In p5| theories with fermions transforming according to higher 
dimensional representations were analyzed. 

In p6| we have put the idea of nonsupersymmetric gauge duality on a firmer ground 
by showing that for certain scalarless gauge theories with a spectrum similar to the one 
of QCD the gauge dual passes a large number of consistency checks. The theory studied 
in ||36ll is QCD with Nf Dirac flavors and one adjoint Weyl fermion. An important feature 
of this theory is that it possesses the same global symmetry of super QCD despite the 
fact that squarks are absent. This means that there are four extra anomaly constraints not 
present in the case of ordinary QCD, moreover we have shown that the potential dual 
can be constructed for any number of colors greater than two. 

The magnetic dual is a new gauge theory featuring magnetic quarks and a Weyl ad- 
joint fermion, new gauge singlet fermions which can be identified as states composite 
of the electric variables, as well as scalar states needed to mediate the interactions be- 
tween the magnetic quarks and the gauge singlet fermions. The new scalars allow for 
a consistent flavor decoupling which was an important consistency check in the case of 
supersymmetry. In p6l| it was also shown that the candidate dual allows to bound the 
anomalous dimension of the Dirac fermion mass operator to be less than one in the con- 
formal window, and there we also estimated the critical number of flavors below which 
large distance conformality is lost in the electric variables. 

In these previous studies we implicitly assumed that there exist nontrivial infrared 
fixed points in the magnetic description matching the electric ones. However, given 
that we do not start in the ultraviolet with a supersymmetric spectrum nor we impose 
supersymmetric coupling relations it is legitimate to ask whether these fixed point exist 
at all in the magnetic description, and if they do what are their universal properties. 

The goal of this work is exactly to establish the existence of these fixed points as- 
suming the most general coupling structure. We will investigate the regime where the 
computations are trustable, i.e. when the magnetic theory features fixed points accessi- 
ble in perturbation theory. We discover, for the first time, a large number of calculable 
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fixed points and classify them by analyzing their basin of attraction. An extremely inter- 
esting result, according to us, is that among these very many fixed points we discover 
supersymmetric fixed points to emerge in the infrared along nonsupersymmetric flows. 

Our results have an impact on the (non)supersymmetric gauge-gauge duality and even 
recent ideas on why the SM can be seen as the magnetic dual description of an electric 
strongly coupled theory featuring no scalars Il37ll . The consequences of this interpretation 
of the SM leads to a novel way to tackle the hierarchy problem and even to shed light on 
the mystery of the observed number of matter generations [37J. 

The phenomenological interest of our studies relies also on the fact that theories similar 
to the ones investigated here and featuring infrared fixed points have been used to 
construct sensible extensions of the standard model of particle interactions of technicolor 
type passing precision data and known as Minimal Walking Technicolor models ||3l|8l. 
Duality has also been used to infer relevant insights on the constraints on these models 
from precision electro weak tests ||38] - |40| . 

II. MAGNETIC SETUP OF QCD WITH ONE ADJOINT FERMION. 

In p6l] we investigated the possible existence of a magnetic dual for the electric theory 
constituted by a scalarless SU{N) gauge theory with Nf Dirac fermions and N larger than 
two, as in QCD, but with an extra Weyl fermion transforming according to the adjoint 
representation of the gauge group. The quantum global symmetry of the electric theory 
is therefore: 

SUdNf) X SURiNf) X Uv{l) X Uaf{1) ■ (1) 

At the classical level there is one more Ua(1) symmetry destroyed by quantum corrections 
due to the Adler-Bell-Jackiw anomaly. Of the three independent U{1) symmetries only 
two survive, a vector like LZy(l) and an axial-like anomaly free (AF) one indicated with 
UafO')- The spectrum of the theory and the global transformations are summarized in 
table H 

We then postulated the existence of a magnetic dual featuring the minimal spectrum 
of composite states and gauge group structure needed to: 

• Classify all of the 't Hooft anomaly conditions and show how to match them, for 
the first time, for any number of colors and flavors; 
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Fields 



A 
Q 
Q 
G„ 



[SU{N)] 



Adj 

□ 
□ 

Adj 



SUiiNf) SURiNf) Uv{l) Uaf{1) 



1 


1 





1 


□ 


1 


1 


N 
N,- 


1 


□ 


-1 


N 
Nf 


1 


1 









TABLE I: Field content of the electric theory and field transformation properties. The 
squared brackets around SU{N) indicate that this is the gauge group. 

• Allow for consistent flavor decoupling both in the electric and in the magnetic 
theory; 

• Ensure duality involution. Meaning that if we dualize once more the magnetic theory 
one recovers the gauge structure of the electric theory. 

The proposed nonsupersymmetric magnetic gauge theory II36II is summarized in Table 
|IIj We have also shown that one has to have X = N f - N. Differently from the dual of 



Fields 



M 



G, 



[SU{X)] 



Adj 

□ 
□ 
1 



□ 
□ 

Adj 



SUtiNf) SURiNf) Uvil) UAFil) 



1 

□ 
1 

□ 



1 

1 

□ 

□ 





Nf-X 



Nf-X 



Nf 



Nf 



Nf-2X 



1 

□ 
1 



□ 
1 
1 



Nf-X 



Nf-X 



X 





Nf-X 
Nf-X 



TABLE II: Field content of the magnetic theory and field transformation properties. 
The four upper fields are Weyl spinors in the (1/2, 0) representation of the Lorentz 
group. The two (/)-fields are complex scalars and G^, are the gauge bosons. 



super QCD we did not impose a supersymmetric spectrum or supersymmetric coupling 
relations. The spectrum above is nonsupersymmetric as it is clear from the fact that there 
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is no complex scalar partner of M. Moreover we showed that we could build the gauge 
singlet states using electric variables which do not contain squarks. Subsequently in [371 
it was argued that one could also add the complex scalar H (appearing now in Table [III]). 
This is possible since it does not affect the anomaly conditions and can be built naturally 
out of the electric fermionic variable as follows [|37|: 

H-QAAQ. (2) 

It is crucial to be able to construct all these states directly from the electric fermionic 
variables. This demonstrates that supersymmetry is not a fundamental ingredient in 
order to construct these states. H in | |37| plays the phenomenologically relevant role of the 
SM-like Higgs, elementary in terms of the magnetic variables. Although when adding 
the new complex scalar field the spectrum of the dual theory looks supersymmetric, 
the full theory is not since the couplings, in the ultraviolet, are not taken to respect 
supersymmetric relations. Assuming the duality to exist we were able to make a number 
of relevant predictions for the nonperturbative gauge dynamics of the strongly coupled 
QCD gauge theory with one fermion in the adjoint representation. The most relevant 
ones being: 

• The anomalous dimension of the mass of the electric f ermions at the lower boundary 
of the conformal window cannot exceed unity. 

• Estimate the size of the conformal window, i.e. the critical number of flavors below 
which the theory looses large distance conformality. 

• When the magnetic theory is used as possible dual of a minimal extension of the 
SM one can argue that this extension can be re-constructed in terms of a strongly 
coupled electric theory featuring only fermionic matter, and furthermore argue that 
the mathematical consistency of the electric dual requires, at least, three generations 
of ordinary matter [|37||. 

What was still missing is the actual existence of the fixed points in the magnetic 
theory. This is what this paper aims to accomplish using perturbation theory. This regime 
is achieved by choosing the number of flavors Nf, number of colors X, and Yukawa 
couplings where perturbation theory holds. 
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Fields 



Am 
1 

M 
H 



[Sii(X)] 



Adj 

□ 
□ 
1 
1 

□ 

□ 

Adj 



SUiiNf) SURiNf) Uv{l) Uaf{1) 



1 
□ 
1 
□ 
□ 
1 
□ 
1 



1 
1 

□ 
□ 
□ 
□ 
1 
1 



Nf-X 
Nf-X 



_X_ 
Nf 

_X_ 

Nf 

Nf-2X 



^ 

N/-X Nf-X 



X 
Nf-X 



r 



Nf 
Nf-X 



X 





Nf 







TABLE III: Field content of the magnetic theory with the addition of the Higgs-field, H 
and the field transformation properties. The four upper fields are Weyl spinors in the 
(1/2, 0) representation of the Lorentz group. 

III. MAGNETIC POTENTIAL OF THE THEORY AND BETA FUNCTIONS 

We start with listing the Yukawa operators: 

£y = yA(p*/^mq + yj(p*Amq + yM(pMq + yM(pMq + yuqHq + h.c. (3) 

(^^-interactions do not affect our results to the perturbative order we are considering and 
therefore we will not include them here. 

This Yukawa-sector contributes to the two-loop beta function of the gauge coupling in 
the following manner: 



dlnfi '^1471)2 '^'(471)4 '^^(471)4 

^ = y C2(G) -^Y, T{r)Nf{r) - ^ J]^ T{s)Nf{s) 

r s 

Pi = ^C^iGf - [^C2(G) + 2C2(r)] T{r)Nf{r) - [^2(0) + 2C2(r)] T{s)Nf{s) (6) 



(4) 
(5) 



^^=^LTr[C2(r)ry/] 



(7) 



d{G) ^ 

where r denotes the representation of fermions and s denotes the representation of the 
real scalars. T(-) is the trace normalization of the group generators, Czi-) is the quadratic 



Casimir of these and d{G) is the dimension of the gauge group. We refer to Appendix 
A for a careful derivation of the Yukawa contribution jSy to the running of the gauge 
coupling. The result is: 



|So + agjSi + («A + a-^ N; + y ^ + anW^ 



(8) 



with 



|6o = 3X - N/, |6i = _ 7^^^ ^ 



3Nf 



X ' 



(9) 



and where we used the notation 



(47t)2 



We now consider the running of the Yukawa couplings. The one loop beta function for 
the Yukawa couplings is given by ||42ll43l : 

(47i)2|6(Y0 = \ \(\(j)Y^ + Y^Y2(r)] + 2Y'^Y+Y'^ + ^ Y'^Tr [y+Y^' + Y]Y^\ - ?>i{C^(r), Y% (10) 

Here Y^ is the Yukawa coupling matrix defined by the particular interaction: 



where roman indices contract over the scalar gauge-flavor overall index and the greek 
indices a, jS are again gauge-flavor indices but reserved for the Weyl fermions ^\) and x- 
Y2(r) is the group invariant: 



Y2(r) = Y]YK 



For each scalar contraction we multiply, row by column, the Yukawa matrices over the 
fermion indices. We again report the derivation of the beta function for each Yukawa 
coupling in Appendix |Aj To the second order in the couplings the set of beta function 



equations reads: 



]6(«a) = lot} 



3X^ — 1 

|6o + «giSi + («A + (X-f) iv/ 



(11) 



3(X2-1) + \ ^ 3X2-1 
— — ax + (ttA + + y ^ + «Mj - 



4N 



f- 



(47Z)2 



j6(aM) = 2aM 



3N/ + X 



X 



2 



x^-iyMy^yAyA 



2X 



— +aA -3a;g 



(12) 



IX (47l)2 



«A + <2!T 



2X 



6a „ C2(n) + X^H 



2X 



/3(tt^) = (yM ^ yM' yi ^ y/i) , 



where: 



iS(«0 = 



2y,- 

(471)^ 



r^(yO • 



(13) 
(14) 

(15) 
(16) 

(17) 



We impose CP invariance and therefore set all the Yukawa phases to zero. We discuss the 



full set of equations, including the phases, in the Appendix A 



In the supersymmetric limit ax = aj = 2a g, aM = a^ = a^ and the beta function system 
simplifies to: 



p{ag) = -2al [|So + ^[(Xg + INjau] 
^{ax) = -2ax [jSo«g] 
|6(aM) = 2«M [{2Nf + X)«M - 4«^^C2(n)] 
|6(«h) = 2«H [{2Nf + X)aH - 4«^C2(n)] , 



(18) 
(19) 
(20) 
(21) 



where ^[ = + {3X^ - 1)N//X = 6X2 - 4JV/X + 2JV//X is the well-known two-loop 
term of pure super QCD, i.e. without the gauge-singlet chiral superfield transforming as 
bilinear with respect to the non-abelian global symmetries and, if present, leading to the 
introduction of a superpotential term. The theory with this chiral superfield included is 
Seiberg's magnetic dual. 

This result shows, as expected, that supersymmetry stays unbroken along the renor- 
malization flow. Also note how the running of the gaugino Yukawa coupling has collapsed 



to the one-loop result for the running of the gauge coupling. These results are in agree- 
ment with the known result from supersymmetry. The all-orders supersymmetric beta 
functions for Seiberg's magnetic dual are: 

^'{au) = au [ym + 2yo] , (22) 

where the one-loop expressions for the anomalous dimensions of the chiral superfields 
read: 

7o = -4C2(n)ag + INfUu + ■■■ 

Ym = 2XaM + • • - (23) 

by which it is readily seen that the one-loop expansion of jS^(«g) and ^^(aM) are in agreement 
with the expressions in equations (|T8|-(21 ). 



IV. MAGNETIC FIXED POINT ANALYSIS 

Just below the critical number of flavors where the theory looses asymptotic freedom 
(j6o = 0) a Banks-Zaks perturbatively stable infrared fixed point (IRFP) emerges once the 
Yukawa interactions are set to zero. The well known expression for the value of the gauge 
coupling at this fixed point is: 

-r = -| ■ p4) 

In this phase the theory shows large distance conformality. We now investigate whether 
the perturbatively stable IRFP persists when we turn on the Yukawa interactions. From 



Eq.([TT]), the perturbative fixed point value for the gauge coupling reads: 

4X^0 + + al){3X' - l)Nf + 2X{a*^ + a*- + 2a*^)N] 
= • ^^^^ 

Perturbative consistency is ensured by the smallness of the Yukawa couplings at the fixed 
point. There is, clearly, a symmetry between the tilded and untilded Yukawa coupling 
constants at the beta functions level. This symmetry persists at the fixed point effectively 
reducing the space of solutions. 
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The fixed points associated to zeros of (iian) are: 

12C2(n)a* - C2(n)(a; + a*-) - Nf{a*^ + a*^ 



0*^ = 



A 



«H = 



2{Nf + X) ■ ^^^^ 

A perturbative nontrivial infrared fixed point for all the couplings is achieved when 
the number of flavors are such that asymptotic freedom is almost lost. To express this 
point mathematically we introduce the small parameter e as follows: 

^o>O^Nf< 3X, Nf = 3X(1 - e) . 

We report in Table IV and Table |V| the explicit solutions for the fixed points as a function 



of e. There are 18 physical solutions for the system ( 11 1-( 16 1, where 10 of these correspond 
to the class a*^ = and the remaining 8 to a*^ 0. The solutions obtained exchanging 
tilded couplings with untilted ones 



a;{ <r^ aj (x 



M 



are indicated in the Tables IV and[V|by an asterisk. Solution (5) in Table [TV| corresponds 
to the perturbative infrared fixed point for super QCD while the last solution (6) in Table 
[Vjcorresponds to the perturbative fixed point for Seiberg's magnetic dual and it is the only 
one with all nonvanishing couplings. Interestingly, when setting to zero = (solution 



7 in Table IV I we find yet another fixed point corresponding to the infrared stable fixed 
point of the non-supersymmetric magnetic dual theory proposed in [|36| and shown in 
Table |llj This lends further support to the proposed non-supersymmetric duality. 



V. STABILITY ANALYSIS 



The stability of the fixed points reported in Tables IV and |V] is determined by first 
evaluating the following matrix 

"■'■=4 (27) 

at the fixed point of interest denoted in brief by g*. If it has only real and positive 
eigenvalues the fixed point is said to be stable. In the Tables [Vl| and [VIl] we present the 
eigenvalues corresponding to the fixed points of Tables IV and |V] respectively for X = 4 
and e — 0.05. For convenience, we have marked the cells with negative eigenvalues in 
gray 
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TABLE IV: Fixed point solutions on the critical surface a*^ = 0. Solution 2,3 and 6 are 

doubled by symmetry property (*). Note that to linear order in e solution 5 is the super 
QCD fixed point. 



(9(e)-expansion 



a, 



_ {3X-Nf)X 
~ N/(7X2-3)-6X3 



a* - - = 



- Xe 
5X2-3 



2W 



2(3X-N/)(X(N/+3X)-3)X 
2N2x(7X2-3)+3N^(5-14X2+X'')-36X3(X2-1) 

12(3X-N/)(3X2-1)X 
2N2x(7X2-3)+3N/(5-14X2+X*)-36X3(X2-1) 



al - = 



_ 2X(2X2-l)e 

~ 5-16X2+11X4 

_ 4X(3X2-l)e 
~ 5-16X2+11X4 



3W 



2(3X-N/)(3N/+X)X 
3N2(13X2-5)-12X4-2N/X(3+11X2) 
, _ 6(3X-N|)(X2-1)X 

~ 3N2(13X2-5)-12X''-2N/X(3+11X2) 

a* - - = 



_ 20Xe 



91X2-51 

6(X2-l)e 
X(91X2-51) 



(3X-Ny)(3N^+2X)X 
2(N2(9X2-3)-6X4-N/X(3+2X2)) 
» ^ 3(3X-Ny)(X2-l) 

2(N2(9X2-3)-6X4-N/X(3+2X2)) 



a - 



= 



llXe 



2(23X2-12) 

_ 3(X2-l)e 
~ 2X(23X2-12) 



6W 



^8 = 



(3X-N/)(3X2+2N/X-3)X 



2(9X3-9X5+n2x(7X2-3)+N^(3-6X2-9X4)) 

^ _ 3(3X-N/)(3X2-1)X 
^ ~ 2(9X3-9X5+n2x(7X2-3)+N/(3-6X2-9X4)) 



a - 



= 



2(3X-N/)(3(N2-1)X-8N/+9N/X2+3X^)X 



- Xe 



2(X2-1) 



a, 



3N3x(13X2-5)-36X4(X2-1)+3N/X(5+18X2-31X4)+2N2(18-59X2+9X4)) 

12(3X-N^)(3X3-2Ny-X+8N/X2) 
3N3x(13X2-5)-36X''(X2-l)+3NfX(5+18X2-31X4)+2N2(18-59X2+9X'')) 
, _ 6(3X-Ny)(N/X-2)(X2-l) 

~ 3N3X(13X2-5)-36X4(X2-1)+3N/X(5+18X2-31X4)+2N2(18-59X2+9X4)) 



_ 2X(19X2-9)e 



41-141X2+100X4 

_ 4X(27X2-7))c 
~ 41-141X2+100X4 

_ (6X2 -4)e 



100X3 -41X 



a- 



KM 



= 



(3X-N/)(3(N2-l)X-4Ny^+6Ny^X2+3X») 
2(9X3-9X5+n3(9X2-3)+N2(X-15X3)+N^(3+6X2-21X4)) 

3(3X-N^)(5NfX2-2Ny-X+3X3) 



«A = 



a- - 



a 



M 



\ 9X3-9X5 +n3(9X2-3)+N2(X-15X3)+N/(3+6X2-21X4) 
^ _ 3(3X-N/)(N/X+3X2-2)(X2-1) 
'^M ~ 2X(9X3-9X5+n3(9X2-3)+n2(X-15X3)+N/(3+6X2-21X4)) 



_ X(16X2-5)e 
~ 2(4X4-5X2+1) 

_ X(18X2-7))e 
~ 4X4-5X2+1 

_ (3X2-l)e 
~ X(4X2-l) 
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TABLE V: Fixed point solutions on the critical surface a*^ 0. Solutions 2 and 3 are 
doubled by symmetry property (*). Note that to linear order in e solution 6 is the fixed 
point of Seiberg's magnetic dual. 



(*) ax <-> a-j^, au ^ oi 



M 



(3(e)-expansion 



_ {3X-Nf){Nf+X) 



4N2X-6X3+N/(X2-3) 



a 



_ (3X-Nf){X^-l) 
H X(4N2X-6X3+N/(X2-3)) 







_ 4Xe 
11X2-3 

_ 3(X2-l)e 
~ X(llX2-3) 



2(*) 



3W 



(3X-N/)(6X(X2-l)+N/(7X2+2XN/-5)) 
2XN2(5X2+4XN/-13)-36X3(X2-1)-3N/(4X2+9X''-5) 
3(3X-N/)(x2-l)(3X2+2XN/-5) 



H X(2XN2(5X2+4XNy-13)-36X3(X2-l)-3N/(4X2+9X4-5)) 



a, = 



12X(3X-N/)(3X2+2XNf-l) 



A 2XN2(5X2+4XNy:-13)-36X3(X2-l)-3N^(4X2+9X4-5) 



a- - = 
AMM 

al = 



(3X-N/)(2X2+8XN/+5N2) 



g 20N3x-12X4+N2(17X2-15)-2N/X{17X2+3) 
3(3X-N/)(2X+5Ny^)(X2-l) 



a' = 



H X(20N3X-12X*+n2(17X2-15)-2N/X(17X2+3)) 
6(3X-N/)(X2-1) 



a', = 



M 20N3x-12X4+n2(17X2-15)-2N/X(17X2+3) 



a* - - = 



_ X(15X2-7)e 
~ 5-26X2+21X4 

_ (9X2-5)e 
~ X(21X2-5) 

_ 4X(9X2-l)e 
~ 5-26X2+21X4 



_ 71Xe 



193X2-51 

_ 51(X2-l)c 
~ X(193X2-51) 

_ 6(X2-l)e 
~ X(193X2-51) 



(3X-N/)(2N/+X)(N/+2X) 
2(4N3x-6X4+n2(7X2-3)-N/X(3+8X2)) 

3(3X-N/)(N/+X)(X2-1) 
H ~ X(4N3X-6X*+N2(7X2-3)-N/X(3+8X2)) 



al = 



a 



a,, = a~ = 



3(3X-N/)(X2-1) 



M 2(4N3x-6X*+n2(7X2-3)-N/X(3+8X2)) 



a* - = 

A,A 



_ 35Xe 
2(47X2-12)) 

_ 12(X2-l)e 
~ X(47X2-12) 

_ 3(X2-l)e 
~ 2X(47X2-12) 



a„ 



{3X-N f){3X{X--l)+2N f{2X-+XNj-l)) 
2(4N/X+3X2-3)(N2X-3X3-N/(X2+1)) 
» _ 3(3X-N/)(N/X-1)(X2-1) 

~ X(4N/X+3X2-3)(N2x-3X3-N/(X2+1) 
3X(3X-N^)(3X2+2XN^-1) 



n = n = 

A X (4N/X+3X2-3)(N2X-3X3-N/(X2+1)) 



a 



MM 



= 



_ X(llX--3)c 
~ 2(5X2-1)(X2-1) 

_ (3X2-l)e 
~ X(5X2-1) 

_ X(9X2-l)e 
~ (5X2-1)(X2-1) 



a„ 



(3X-N/)(2N/+X)(N2+2N/X+3X2-3) 



2(4N*X+9X3-9X5+n3(X2-3)-N2x(3+23X2)+3N/(1+4X2-9X*)) 
3(3X-N/)(N/X-N2+1)(x2-1) 

X(4N*X+9X3-9X5+n3(X2-3)-N2x(3+23X2)+3N/(1+4X2-9X*)) 



a' 



3(3X-Nf){X+2Nf){3X^+XNf-l) 



X 4N4X+9X3-9X5+JV3(X2- 



3) 



-N2x(3+23X2)+3N/(1+4X2- 



9X4)) 



a~ = 



3(3X-N/)(3N/X+3X2-2)(X2-1) 



7Xe 
2(X2-1) 

X 



7Xe 
X2-1 



lal + Oie^) 



M 2X(4N|X+9X3-9X5+n3(X2-3)-N2x(3+23X2)+3N/(1+4X2-9X4)) 
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i=a*^ + Oie^) 



TABLE VI: Eigenvalues corresponding to the fixed points in Table IV for X = 4 and 
e = 0.05 case. 



# 




V2 


V3 


V4 


V5 


V6 


1 
i 








u.uio 


U.UiO 




9 
Z 


u.zu 






U.UID 


U.UUDO 


n nni i 

U.UUi i 


a 

J 














4 


-0.064 


0.038 


0.030 


-0.027 


-0.023 


0.0039 


5 


0.30 


0.092 


0.063 


-0.060 


-0.058 


0.0089 


6 


0.20 


-0.062 


-0.043 


0.031 


0.021 


0.0065 


7 


3.0 


1.8 


1.4 


1.1 


-0.14 


0.033 



TABLE VII: Eigenvalues corresponding to the fixed points in Table V for X = 4 and 
e = 0.05 case. 



# 




V2 


V3 


V4 


V5 


V6 


1 


0.12 


-0.067 


-0.067 


-0.0073 


-0.0073 


0.0058 


2 


0.37 


0.17 


-0.078 


0.026 


0.012 


0.00038 


3 


0.12 


-0.067 


-0.064 


0.015 


-0.0075 


0.0058 


4 


0.12 


-0.074 


-0.053 


0.015 


0.014 


0.058 


5 


0.73 


0.26 


0.18 


0.16 


-0.081 


0.018 


6 


21 


5.4 


4.0 


4.0 


0.77 


0.052 



There is a rich structure of fixed points classified according to the relevant directions 
clearly visible from the list of associated eigenvectors reported in the Appendix |Bj The 
dimensions of the critical surfaces are dictated by the number of positive eigenvalues. 
Once the unstable directions, in the coupling space, are removed the remaining subset of 
couplings generate a critical rf-dimensional surface with d given by the number of positive 
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eigenvalues containing the nontrivial fixed points. As an example consider the first fixed 



point of Table VI Here we have a critical line with an infrared stable fixed point and five 
unstable directions. The critical surface in this case is just a line. 



Another interesting example is constituted by the fixed point solution 7 of Table IV 



In this case we find a five-dimensional critical surface and one unstable sixth direction 
parallel to the i/h axis. We plot in Figure [T] the projections of the renormalization group 
flow around this solution for a subset of two couplings at the time. The other filled small 
circles in the figure correspond to the projections on the chosen planes of fixed points 
such as the Seiberg's magnetic dual, indicated by the SUSY label on the plot, or the super 
QCD Banks-Zaks fixed point labelled by susyBZ. 

The fixed point featuring all positive eigenvalues is the one corresponding to Seiberg's 
magnetic dual. 



VI. PHYSICAL RESULTS AND CONCLUSIONS 



We uncovered the full spectrum of perturbative fixed points associated to an SU{X) 
nonsupersymmetric gauge theory featuring scalars and fermions, of the type summarized 



in Table III Although the nontnteracting field theory has a supersymmetric looking spec- 
trum we did not assume the bare couplings to respect supersymmetry. We then analyzed 
the fixed points by studying the beta functions to second order in the couplings. To this 
order, which is renormalization scheme independent, we discovered supersymmetric and 
nonsupersymmetric fixed points. We further analyzed their stability and discovered that: 

• When all bare couplings are nonzero the theory flows to the Seiberg's magnetic dual 
fixed point. This occurs on supersymmetric and nonsupersymmetric renormaliza- 
tion flow directions. 

• When the Yukawa couplings i/m/ Vh to zero the theory flows to 
the super QCD fixed point, as above, independently if the renormalization flow is 
supersymmetric. 

• When all the Yukawas are set to zero we achieve the perturbative nonsupersym- 
metric fixed point in the gauge coupling. 
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FIG. 1: Projections of RG flows for X = 4 and e = 0.05 in the planes (from left to right 
and top to bottom): (i/m, i/h), ig, Vh), {y\,yH) and (i/m/ Vm)- ^^her couplings are kept 



fixed at the fixed point value of solution 7 given in Table IV Besides solution 7 clearly 
indicated in the plot we also plot the projections of the Seiberg's magnetic dual 
indicated by the SUSY label on the plot, and the super QCD perturbative fixed point 
indicated by susyBZ. 



• When xfM and are nonzero and i/h = we discover a new nonsupersymmetric 
fixed point. This can be identified with the magnetic dual fixed point for nonsuper- 
symmetric gauge theories proposed in p6|. Another relevant fixed point emerges 
when = = and yu + 0. 
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• All the other fixed points emerge similarly by setting to zero different Yukawas. 
The generic conclusions we draw from our results are: 

• Supersymmetry naturally emerges as a fixed point theory from a nonsupersymmet- 
ric Lagrangian, meaning without fine-tuning the bare couplings. This relevant result 
demonstrates that supersymmetry can be viewed as an emergent phenomenon in 
field theory. Our results, among other things can be used to argue that one does not 
need to fine-tune the bare couplings when performing Lattice simulations aimed to 
study supersymmetry on the Lattice ||44ll45l . 

• Duality is not a prerogative of supersymmetry given that the theory features many 
nonsupersymmetric fixed points, depending on the details of the interactions. 
Therefore these new fixed points can describe the nonperturbative physics of electric 
dual gauge theories like the one envisioned in Il36|| . 

The potential physical applications of our results range from (supersymmetric) exten- 
sions of the SM to providing new ideas on how to investigate nonperturbative supersym- 
metric dynamics via first principle lattice simulations. 
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Appendix A: Derivation of the beta functions 

We follow the notation of Machacek and Vaughn [|4Tl |42| . who derived the general 
expressions to two-loop order for the running of the gauge coupling and Yukawa cou- 
plings for a general gauge theory with real scalars and Majorana fermions. In Ref. [43l 
is showed how to use the expressions of Machacek and Vaughn for complex scalars and 
Weyl fermions. We need only to do a slight rewriting of Xy given in Eq. ^ to follow their 
notation: 

£y = yA(p*rqA^^ - yjqTYA^^ + yM(pMq + V^Mq + yHHqq + h.c. 

= yj^cp^^qaKr + ykar;^"^ + U^cPjMq" + UJ^'Mq, + V^Hq^f + h.c, (Al) 

where a = 1,. . . , d{G) is the gauge index reserved for the adjoint Majorana fermion with 
d{G) the dimension of its representation, greek gauge indices «, j6, . . . are reserved for 
the Weyl fermions and roman gauge indices i, j, . . . are reserved for the complex scalars. 
Flavor indices have been suppressed. The Yukawa matrices are defines as follows, once 
we take flavor indices into account with 1,1' denoting the indices of SlI{Nf)i and r,r' 
indices of the SU{Nf)^ global symmetries: 





(A2a) 


W = -y-(T") ^(5'' 

^ ar,a i/A^ 'a '-'r ' 


(A2b) 




(A2c) 


= m^pA" 


(A2d) 




(A2e) 



where T" are the generators of the gauge group. The minus sign of y^-term has been 
chosen to facilitate the connection with supersymmetry, i.e. for SUSY yA = yj = i V2g. 

We will first derive the Yukawa contribution to the 2-loop beta function, given in eq. 
Q-Q, in particular: 

where C2{r) is the quadratic Casimir of the spinor representations r and d{G) is the dimen- 
sion of the gauge group. It is instructive to consider the contribution from the Yukawa 
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sector to the two-loop beta function of the gauge coupling diagrammatically. There is 
only one type of diagram that has non-vanishing contribution, which is: 




where the solid line represents spinor fields, and the dashed line represents the scalar(s) 
coupled via the Yukawa interactions. Note that the scalar-gauge interactions do not 
contribute. 

The fermion M is not coupled to the gauge fields. Thus we get: 

+ C2(^ [VY^j + WU,^ + VV^]\ 
= (^A + aj)T{n)Nf [C2(n) + €2(0)] + {au + «m + 2«h)N- 



-2C2(D)X 

f d{G) 



N/ 3X2-1 /«M + «M \^,2 /A.N 



where we used the notation 

_ \y,\' 

(471)2- 

We now consider the running of the Yukawa couplings. The one loop beta function for 
the Yukawa couplings is given in | |42ll43ll : 

(47i)2j3(y0 = ^ [Yl{r)Y> + YiY2{r)] + 2Y''Y^.Y'' + ^Y'^Tr [yIY' + Y+Y*^] - 3g^{C2{r), Y'}. (A5) 

where Y2(r) is the group invariant: 

Y2(r) = YlY' (A6) 

It is again instructive to consider each term in the beta function in terms of the diagrams 
they originate from. The first two terms in the bracket corresponds to the fermion leg 
self-energy contribution from the Yukawa sector, and the expression tells us that they 
each contribute with a factor of | to the beta function: 
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The next term in the beta function is the vertex correction from the Yukawa sector: 




The fourth term is the scalar leg self energy contribution from the Yukawa sector: 




Note the different directions of the arrow in the fermion loops. For Dirac fermions the 
above contribution comes naturally with a factor of two. Finally the gauge sector has a 
non-vanishing contribution only to the fermion leg self energy, and each term contributes 
with a factor of -3: 
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We thus jfind the one-loop coefficient of the beta function for Y' to be: 

+ y'^Tr [y^r] + y'^Tr [W'U^] - 3g^ {C2{q)Y' + y'C2{A)} . (A7) 

Inserting the definitions of the Yukawa-matrices, we find: 

i^nflHYi) = {C2(n) I yA P +Nf \ yu ? +Nf \ yn f +T(n)N/(| y,f + \ yj |^)) 

+C2(n) I yA \' +Nf I yu I' -3f ma) + C2(G)}] ^ - lNfy'^y^yj(r)Jd\ 

Note that (r")^^6|' in the last term is the matrix part of Y^, which cancels on both sides of 
the equation once we do the rewriting j3(y^) = {T'^)Jb^^^{y\). 
The beta function for Y^ is obtained via: 



yA <^ yA yM <^ yM 



^{W) and j3(L[^) are related by a similar transformation. We compute the former: 

(47i)2i3(LJ0 =^ [{U\{M) + U\{M)] W + W {Ua® + yz® + ^2®)] + 2ry;LJ' 
+ U'^Tr [UIU>] + U'^Tt [YlYi] - 3g^WC2(^. 
Using the expression for the matrices we deduce: 



mo = 



+NfaM + C2{a)a, - 3a,C2{a)] W - ^^yiy^yMSisf 6^', (A9) 

where we used ai =| y, p /{4cn)^. As before is the matrix part of W, which cancels 

on both sides of the equation. 

The renormalization of the yn coupling reads: 



(A8) 



{4nfm = 



+ 



■ [Yliq) + Uliq) + Vliq)] V + v[U2(q) + Y2(^ + V2(q)} 
VlT [V'V] - Sg" {C2{q)V + VC2m 



^ N/ + y — 6a, 



)c2(n) + X«H 



V 



(AlO) 
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We summarize now the results for beta functions, after having cancelled the matrices: 



|So + Ug^i + NfT{n) [C2(n) + C2(G)] {ua + aj) + 



C2in)XNj 



d{G) 



(Alia) 



|6(y,) =1 I e'^'^ 



3^ ^ ^ Tin)Nf 
^C2{n)aA + — - — 



(47l)2 

3 



(^A + «x) + \^-^ + «m] - 3«, {C2(n) + C2(G)} 

(Allb) 



^Nf \aM + y) + f («M + «m) + <^2(n) + ^a] - 3«^C2(n) 



i8(yM) =1 vm I f^''" 

_ 2C2(n)^^4^^^^e'(^M+A0.) 



(47Z)2 



Nf /a,\ + aj 
—{au + «M + 2«h) + C2(n) I — ^ I + X^h 



where we have made explicit the phases of the Yukawa's couplings: 

Vi =1 Vi I e'^'' 

A0M = ~ ^^A = ~ ^A- 



(Allc) 

(Alld) 
(Alle) 



Note also that we define 



(47l)2 



^ 7£^^(y<), 



where y, is the complex conjugate of y,. 

We have considered the case of real couplings in the main text. 

Appendix B: Tables 



We provide here the eigenvectors associated to the eigenvalues given in Tables VI and 



VII Shaded cells indicate unstable directions, i.e. corresponding to negative eigenvalues. 
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TABLE B.l: Eigenvectors associated to v/ given in Table VI and corresponding to the 



fixed points 1,2 and 3 of Table IV 





Vm 


Vm 


1/A 


3/1 


g 


yn 


Vl 











1.0 








V2 








1.0 











V3 

















1.0 


V4 





1.0 














V5 


1.0 

















V6 














1.0 






VTl2 


1/M 


Vm 


1/A 


yj 


g 


yn 


Vl 











-53. 


1.0 





V2 








1.0 











V3 

















1.0 


V4 


1.0 

















V5 











1.8 


1.0 





V6 





1.0 















VTl3 


1/M 


3/m 


yA 


y-A 


g 


yn 


Vl 








1.0 











V2 











1.0 








V3 





-46. 








1.0 





V4 

















1.0 


V5 


1.0 

















V6 





0.61 








1.0 






TABLE B.2: Eigenvectors associated to V; given in Table VI and corresponding to the 



fixed points 4 and 5 of Table IV 



Vll4 


1/M 


J/m 


yA 


yj 


g 


yn 


Vl 








1.0 


1.0 








V2 


-23. 


-23. 








1.0 





V3 


-1.0 


1.0 














V4 

















1.0 


V5 








-1.0 


1.0 








V6 


0.58 


0.58 








1.0 






Vll5 


1/M 


J/m 


1/A 


yj 


g 


yn 


Vl 








-22. 


-22. 


1.0 





V2 








-1.0 


1.0 








V3 


-1.0 


1.0 














V4 


1.0 


1.0 














V5 

















1.0 


V6 








1.5 


1.5 


1.0 
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TABLE B.3: Eigenvectors associated to v, given in Table VI and corresponding to the 



fixed points 6 and 7 of Table IV 



VTl6 


Vm 


Vm 


Va 


yj 


g 


VH 




-2.2 








-49. 


1.0 










0.22 


1.0 











■ 


1 














1.0 




-40. 








13. 


1.0 





V5 





-0.73 


1.0 











V6 


0.47 








1.8 


1.0 






VTl7 


1/M 


Vm 


Va 


3/1 


8 


yn 




-0.035 


-0.035 


-3.1 


-3.1 


1.0 





V2 


-0.41 


0.41 


-1.0 


1.0 








V3 


0.40 


-0.40 


-1.0 


1.0 








V4 


-2.0 


-2.0 


1.7 


1.7 


1.0 





V5 

















1.0 


V6 


0.64 


0.64 


1.6 


1.6 


1.0 






TABLE B.4: Eigenvectors associated to V; given in Table VII and corresponding to the 
fixed points 1 and 2 of Table |Vj 



Villi 


Vm 


Vm 


Va 


yj 


g 


yn 


Vl 














-0.052 


1.0 


V2 











1.0 








V3 








1.0 











V4 





1.0 














V5 


1.0 

















V6 














1.1 


1.0 



VIll2 


1/M 


yu 


Va 


y-A 


g 


yn 


Vl 











5.2 


-0.28 


1.0 


V2 











-1.3 


-0.057 


1.0 


V3 








1.0 











V4 





1.0 














V5 











2.0 


1.2 


1.0 


V6 


1.0 
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TABLE B.5: Eigenvectors associated to v, given in Table VII and corresponding to the 
fixed points 3 and 4 of Table V 



VII|3 


Vm 


Vm 


Va 


yj 


8 


yn 







0.15 








-0.054 


1.0 


V2 








1.0 











V3 











1.0 








V4 





-6.0 








0.094 


1.0 


1- 


1.0 

















n 





0.47 








1.2 


1.0 



VIII4 


Vm 


Vm 


1/A 


yj 


8 






0.16 


0.16 








-0.057 


1.0 


V2 








1.0 


1.0 








V3 








-1.0 


1.0 








V4 


-1.0 


1.0 














V5 


-2.9 


-2.9 








0.093 


1.0 


V6 


0.49 


0.49 








1.2 


1.0 



TABLE B.6: Eigenvectors associated to V; given in Table VII and corresponding to the 
fixed points 5 and 6 of Table |Vj 



VII|5 


Vm 


yu 


1/A 


yj 


g 


yn 










3.9 


3.9 


-0.54 


1.0 


V2 








-0.93 


-0.93 


-0.076 


1.0 


V3 








-1.0 


1.0 








V4 


-1.0 


1.0 














V5 


1.0 


1.0 














V6 








1.8 


1.8 


1.3 


1.0 



VII|6 






1/A 


yj 


g 


yn 


Vl 


0.41 


0.41 


5.1 


5.1 


-7A 


1.0 


V2 


-0.44 


0.44 


-1.0 


1.0 








V3 


0.37 


-0.37 


-1.0 


1.0 








V4 


1.7 


1.7 


-3.8 


-3.8 


-1.2 


1.0 


V5 


-0.38 


-0.38 


-0.26 


-0.26 


0.072 


1.0 


V6 


1.1 


1.1 


2.8 


2.8 


1.9 


1.0 
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